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Abstract 

The weak form factors of the nucleon, including the induced pseudoscalar 
form factor and second class terms, are constrained using a microscopic 
calculation of the weak capture process 3 He(/i~, u^B.. The calculation is 
parameter free, and yields a rate of 1499(16) Hz, in agreement with the re- 
markable experimental measurement 1496(4) Hz. The nuclear wave functions 
are obtained using the EIHH method with the Argonne Vis nucleon-nucleon 
potential and the Urbana-IX three nucleon force. The weak currents in the 
nuclei are described using HB^PT formalism. The induced pseudoscalar form 
factor is found to agree with HB^PT prediction. The result is compatible 
with vanishing second class currents, with the tightest constraint to date on 
the conservation of vector current (CVC) hypothesis. 
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1. Introduction 

The weak process in which a muon is captured by a nucleus provides 
an experimental hatch to various aspects of the fundamental forces. For 
decades, it has been used to constrain the properties and symmetries of the 
weak interaction, and to probe the structure of the nucleus at relatively large 
momentum transfer \q\ ~ = 105.6 MeV [lj. 
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This sizable momentum transfer enhances the effect of terms proportional 
to q. One example is the induced-pseudoscalar form factor. This form factor 
has been the target of numerous studies in the past, which have revealed a 
contradiction between experimental and theoretical estimations. These am- 
biguities have induced a new experimental effort by the MuCap collaboration 
to measure the muon capture rate on protons (fi~p) |2|, which have resolved 
the contradiction. By now, the MuCap measurement [2j has already reached 
a ±2.4% determination of the capture rate, and aims to ±1%. This leads to 
the tightest experimental bound on the induced-pseudo scalar form factor of 
the nucleon, constraining its value to ±15% 0]. 

Two additional poorly known form factors are the second class terms, in 
Weinberg's classification [3|, who also assumed their vanishing. The vector 
second class term is also required to vanish by the conservation of vector cur- 
rent (CVC) hypothesis. These terms are demanded by Lorentz covariance, 
but change sign under G-parity transformation. However, this behavior only 
suppresses their value, due to the fact that isospin symmetry is only an 
approximate symmetry of the strong force. The current experimental de- 
termination of both terms has not reached the sensitivity needed to test the 
theoretical estimations, and is still consistent with vanishing form- factors [I| . 

Further reduction of the uncertainties in these three form-factors demands 
sub-percentage experimental accuracy, which is hard to achieve in fi~p pro- 
cess, since its rate is smaller than the free muon decay rate by a factor bigger 
than 500, and since it results in the emission of neutral particles. These obsta- 
cles can be removed for heavier nuclei, as the capture rate is very sensitive to 
the nuclear charge, scaling as Z A . Alas, due to the strong correlation between 
nucleons, theoretical microscopic studies with sub-percentage precision, are 
possible only in very light nuclei. 

However, one capture process fits the theoretical limitations, and has been 
measured to a very high accuracy. The capture of muon on 3 He which results 
in a triton, 

(T ± 3 He -> ^ ± 3 H, (1) 

with a measured capture rate of ± 3 He — > ± 3 H) cxp = 1496(4) Hz, i.e. 
a ±0.3% precision |5j. This precision measurement has already induced a 
number of theoretical works [6 , among them are also microscopic theoretical 
studies of the reaction 0,0,0. The main conclusion of these works is that 
theoretical evaluation should include both state of the art description of the 
nuclear states, and a correct description of the weak interaction between 
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the muon and the nucleus, including the interaction with meson exchange 
currents (MEC) in the nucleus. Though these studies have achieved very high 
accuracy in the nuclear sector, they missed an important, recently discovered, 
ingredient. Lately [9], the electroweak radiative corrections to the capture 
process have been calculated. The results indicate that the enhancement 
factor due to this effect is RC(He) = 0.030(4). The authors of Ref. Q 
demonstrate that this destroys the good agreement of the theoretical study 
of reaction (flT) , done in Ref. [8J. In light of this development, a new theoretical 
evaluation of the process is called for. 

In the current work we cope with this challenge. We study theoretically 
reaction ([T]), and use it to put constraints on the weak form factors of the 
nucleon, including the second class terms, which was not done previously in 
microscopic calculations. We use a hybrid approach, which has been proved 



efficient in describing different weak processes with A = 2 — 4 nuclei [22], [23 
using phenomenological Hamiltonian for the nuclear states, and heavy baryon 
chiral perturbation theory to describe the weak interaction of the muon with 
the nucleus. The latter puts on common ground the single nucleon current 
and meson exchange currents, a fact which increases the reliability of the 
calculation. In addition, the calculation is parameter free, thus can be used 
to make predictions. 

The letter is built as follows. In the next section we outline shortly the 
standard model formalism for calculating the capture process. The solution 
of the nuclear problem is described in Sec. 3, followed by the derivation of 
the weak currents in the nucleus in Sec. 4. In Sec. 5 we give the theoretical 
capture rate and present an error estimation on the results. We discuss the 
consequences of the results on the weak form factors in the last section. 



2. Theoretical Formalism 

We start with a brief reminder of muon capture process, and the formal- 
ism used in the calculation. The muonic atom, a bound state of a muon and a 
nucleus, is unstable. It has two main possible decay schemes, either through 
free muon decay to lighter leptons, or through a muon capture by the nu- 
cleus. The capture is a weak process, where the negative muon interacts with 
the nucleus through the exchange of heavy W~ boson. As the momentum 
transfer in the process is much smaller than the mass of the W~ boson, the 
weak interaction Hamiltonian is given by Hw — — G ^ d ^ / d 3 xj+(x) J -/i (af), 
where G = 1.166371(6) x 10~ n MeV -2 is the Fermi coupling constant (lo| . 



3 



V U d = 0.9738(4) is the CKM matrix element mixing u and d quarks involved 



in the process 



[x) is the lepton charge raising current, and J M is the 



nuclear charge lowering current. 

It is straight forward to evaluate the lepton current. We first note that 
since Za <C 1 (Z is the nuclear charge, and a is the fine structure constant) 
the muon bound in the atom can be regarded as non-relativistic. Moreover, 
the initial (muonic) atom has a Bohr radius much larger than the nucleus 
radius, thus penetration operators are negligible for the needed accuracy 



The discussion in Ref. [13j shows that in these conditions, one can 
approximate the lepton current as a current of point like Dirac particles, 
whose states are described by plane waves, multiplied by a correction factor. 
This factor takes into account the initial bound state of the muon in the atom, 
and the charge distribution of the nucleus. For 3 He, this factor is calculated in 

av\2 _ j^ (ZaM r ) 3 



Ref. [12|,[13j: 



where TZ = 0.979, and M r = (M^+m- 1 ) 



is the reduced mass of the muonic atom. 



We use the Golden rule to write the capture rate [12 . 
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where J3n e = \ is the total angular momentum of the 3 He. The effects of the 
nuclear interaction are embedded in the nuclear matrix element Tn, which 
can be written using multipole decomposition: 
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Cj, Lj, Ej, Mj are the Coulomb, longitudinal, transverse electric and 
transverse magnetic multipole operators of angular momentum J, built from 
the charged nuclear current. One should notice that since the total angular 
momentum of the 3 He and 3 H is |, and both posses positive parity, only 
some multipoles survive: Cq , Lq, Cj 4 , L±, E±, M\ . The superscript A (V) 
stands for operators of axial (vector) symmetry. As the vector current 
satisfies CVC, the vector Coulomb and Longitudinal reduced matrix elements 
are related: < >= < Cj >. 
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Method 


Binding Energy [MeV] 




a He 


EIHH 


8.471(2) 


7.738(2) 


CHH 


8.474 


7.742 


FY 


8.470 


7.738 


Experiment 


8.482 


7.718 



Tabic 1: Binding energies of 3 H and 3 He calculated using AV18+UIX Hamiltonian model 
compared to the same calculation done by using FY equations and CHH method (20| . 
For the EIHH calculation, the number in parenthesis indicates the numerical error. Also 
shown are the experimental values. 

From this discussion, it is clear that the needed information is the struc- 
ture of the weak currents in the nucleus, and the wave functions of the 3 He 
and triton. We will discuss these issues in the following two sections. 



3. The Nuclear Wave Functions 



The evaluation of in Eq. (j2J), demands the solution of the three-body 
nuclear problem, for the ground states of the triton and 3 He. We solve 
the Schrodinger equation microscopically using the effective interaction in 
the hyperspherical harmonics (EIHH) approach 14|, as implemented in the 
nbody fortran codel5|]. In a previous study of reaction (pQ), Marcucci et 
al. js| have shown that the capture is essentially independent of the nuclear 
force, as long as it describes correctly the binding energies of the trinuclei. 
Thus, the nuclear Hamiltonian is taken as the nucleon-nucleon potential Ar- 
gonne vis (AV18) jl6j with the Urbana IX (UIX) [l7| three nucleon force. 
This Hamiltonian has been used successfully to reproduce the spectra of the 
trinuclei as well as othe r lig ht nuclei [17|, and also electro- weak reactions 



with light nuclei [18|, [19|, [26( . Table [T] shows a comparison of our numerical 
results for the binding energies of the trinuclei with the experimental mea- 
surements. Also shown is a comparison to the calculation made using two 
other ab-initio methods, solving the Fadeev-Yakubovski (FY) equations, and 
using the Correlated Hyperspherical Harmonics (CHH) method 2^ . 
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4. Weak Currents in The Nucleus 



The main difference between previous works and the current one, is in 
the details of the nuclear current. The formal structure of the nuclear charge 

lowering current is dictated by the Standard Model: J~ — ^ ^ + J^j , 

where the superscript A (V) stands for current with axial-vector (vector) 
symmetry. r_ is the isospin lowering operator. The axial and vector currents 
are more complicated, as they are affected from the strong interaction, which 
governs the dynamics in the nucleus. To the best of our knowledge the 
fundamental theory of the strong interaction is QCD. Thus, in principle the 
currents should be extracted from the QCD lagrangian. This is, however, 
impossible due to the non- perturbative character of QCD at low energy. 

A possible solution to this problem is found in an effective field theory 
(EFT) approach to QCD [2lj], that is x?T. In x?T, one uses the fact that the 
QCD lagrangian is chirally symmetric in the limit of massless up and down 
quarks, i.e. it is invariant under global SU(2)l x SU(2)r transformations. 
The absence of parity doublets in the low mass hadron spectrum shows that 
the axial symmetry is spontaneously broken at low energies. The pions are 
identified as the Goldstone bosons of the chiral symmetry breaking, and their 
mass is interpreted as a result of the fact that the u and d quarks have mass, 
albeit small. Thus, %PT constructs a low energy lagrangian, that consists 
of nucleons and pions, and posses the symmetries of QCD. Furthermore, 
Weinberg 21| has given a recipe for organizing this lagrangian in terms of 
(Q/\y, where Q is the typical momentum in the process (about 100 MeV 
for muon capture), or the pion mass, A is of the order of the EFT breakdown 
scale, and v > 0. An additional simplification is due to the large nucleon 
mass, which is of the order of the chiral symmetry breaking scale, which 
allows non-relativistic expansion of the lagrangian, the so called heavy baryon 
xPT (HBxPT). The nuclear currents, from this point of view, are Nother 
currents derived from the axial and vector symmetries of this lagrangian. In 
the last two decades, a huge amount of work has been done to derive the 
nuclear currents from the %PT lagrangian. 



The nuclear currents are derived in Ref . |22|, |23( , from a ;\TT lagrangian. 
The currents are expanded to next-to-next-to-next-to-leading order. As 
expected, one finds a nucleonic current, i.e. the impulse approximation, 
and meson exchange currents (MEC). The nucleonic current achieved in this 
formalism is identical in its form to the usual impulse approximation (IA). 
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Its vector part takes the form: 



jV{lk)=u{p>) 
whereas the axial part is 

% (IA) = -u(p') 



Fv(q 2 )r + ^F M tf)a^q u 



u{p), (4) 



u(p). (5) 



Here, is the nucleon mass, is the muon mass, and u(p) is the Dirac 
spinor of the nucleon of momentum p. In order to keep with the power 
counting of HB%PT, we expand Eq. (HHHD in powers of 1/M N , up to 0(M~^ 3 ). 
The currents contain four form factors. Fy and Ga are the vector and axial 
form factors, Fm is the weak magnetism form factor, and Gp is the induced 
pseudo-scalar form factor. In this order of %PT, the form factors contain 
one-pion-loop correction reflected in their q 2 dependence. The first three 
form-factors are very well determined experimentally. The vector and weak 
magnetism are just isospin rotations of the electro-magnetic form factors, 
-FV(O) = 1 and Fm(0) = 3.706, and are extrapolated to the kinematics of 
reaction i.e. to q 2 = -0.954m;;, F v = 0.974(1) and F M = 3.580(3). The 

2 

momentum dependence of the axial form factor is (^(g 2 ) = «7a(1 + if? 2 )^ 
with gA = 1.2695(29) and the axial radius of the nucleon r\ = 0.43(3) fm 2 , 
thus G A (-0.954m 2 ) = 1.245(4). 

For the induced pseudoscalar coupling, HBxPT prediction to one loop 
corrections [24| coincides with the well known Adler-Dothan formula: 

= 2m f^U _ 1 MnT 2 (6) 

m* — q 2 3 

with g npn = 13.05(20) and f v = 92.4(4) MeV. In our case g p = 7.99(20). 

The single nucleon currents are invariant under G-parity transforma- 
tionfl and were thus classified by Weinberg j3j as first class currents. In 
principle, the electro-weak theory does not exclude the possibility of second 
class currents, which change sign under these transformations. It is clear 
that G-parity breaking currents can rise from the fact that isospin symmetry 
is only approximate, i.e. of the order of Int^H^i . Using general symmetry 



The combination of charge conjugation and a rotation in isospin space. 
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arguments, their contribution to the single nucleon currents of Eq. (JUG)]) can 
be written as 5Jj(IA) = ■^-q'*, and 5J^(IA) = -i^-a^q^. In addition, 
the term added to the vector current breaks the well known conservation of 
vector current (CVC) hypothesis. We will test the constraints reaction (pQ) 
can put on these currents in the discussion section of the paper. 

As stressed above, incorporation of MEC is essential for a percentage 
level prediction of the capture rate. In the HB^PT formalism MEC appears 



at v — 2 [22j. In configuration space, the MEC are Fourier transform of 
propagators with a cutoff A. This leads to a cutoff dependence, which is 
renormalized by a cutoff dependent counterterm. Due to the limited scope 



of this letter, we refer the reader to Ref . [22j, [23j , for the explicit form of the 
MEC operators. For the muon capture process, all low-energy coefficients 
in the MEC can be determined from pion-nucleon scattering, except for one 
counterterm d r (A), which characterizes the strength of a two-nucleon contact 
term. 

We fix d r (A) by reproducing the experimental triton half-life of 12.264 ± 
0.018 years, which corresponds to strength of 0.6835 ± 0.001 0, 00. 



Due to the fact that the value used here for the triton half life is used for 
the first time, we use a conservative error estimation due to it, multiplying 
by a factor of 2 the quoted experimental error bar (i.e. we use 0.002). The 
resulting cutoff dependence of d r (A) is: 

d r {A = 500 MeV) = 1.05(6)^0)^ 

d r (A = 600MeV) = 1.82(7) t (l)jv (7) 

4(A = 800MeV) = 3.88(9)^2)^ 

The first error is due to the triton half life, while the second is due to numerics. 

This concludes the nuclear current needed for the calculation, and speci- 
fies the uncertainties in it. 

5. Theoretical Capture Rate and Error Estimation 

The results for the nuclear matrix element of the muon capture process, 
r^, are listed in table [2j The results show a 9% effect due to the MEC 
contribution, and a small effect due to the cutoff dependence of the HB^PT, 



2 This value is slightly different than the one used by Ref. 22|,|26|], consequently changing 
d r (A) calibration. 
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IA 


Total 




A = 500 


A = 600 


A = 800 


Et 

my 

C? 


0.5612 
0.1134 
0.2777 
0.0030 
0.3329 


0.5778 
0.1292 
0.2983 
0.0030 
0.3329 


0.5756 
0.1312 
0.3012 
0.0030 
0.3329 


0.5745 
0.1337 
0.2985 
0.0030 
0.3328 


r^v 


0.6499 


0.7060 


0.7078 


0.7085 



Table 2: The nuclear matrix element Tjv and its different multipole contributions, for the 
impulse approximation (IA) and the calculations which include MEC (Total) for different 
cutoff values A (in units of MeV). 

within 0.3%. One can average this dependence, and arrive at the prediction: 
Tn = 0.7075(10), when using the nominal values of the parameters through- 
out the letter. It is worthwhile noting that the relative contribution of the 
MEC to this process is almost three times bigger than the MEC contribution 
to the triton half life. The extremely weak cutoff dependence shows that the 
essential physics is captured in the HB%PT operators. 
Thus, our final prediction for the capture rate is 

T = 1499(2) A (3) NM (5) t (6) RC Hz, (8) 

where the first error is due to the HB^PT cutoff, the second is due to uncer- 
tainties in the extrapolation of the form factors to finite momentum transfer, 
and in the choice of the specific nuclear model, the third error is related to 
the uncertainty in the triton half life, as reflected in Eq. ([7]), and the last er- 
ror is due to theoretical uncertainty in the electroweak radiative corrections 
calculated for nuclei 0|. This sums to a total error estimate of about 1%. 

The error estimation due to the choice of the specific nuclear model was 
discussed by Marcucci et. al.[8|. By considering different force models, they 
found that using different nuclear potentials does not have a substantial effect 
on the capture rate, as long as the calculation reproduces the binding energies 
of the trinuclei. Their estimate for the uncertainty resulting from this was 
about 2 Hz. As the evaluation of HB^PT based potentials evolves 29], one 
would be able to use a nuclear model of the same microscopic origin as the 
currents, which could merge the nuclear model error and the A cutoff, and 
possibly reduce the estimated uncertainty. 
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The experimental error in the triton half life, as mentioned earlier, is 
multiplied by 2, which will be reduced in the future, when the current mea- 
surement will become standard. 

The radiative corrections, which are the source of the largest contribution 
to the error estimation, are not taken into account in previous studies jsl, 0]. 
One could reduce the large uncertainty in this contribution, by extending 
the work in Ref. [sj to include higher order effects and incorporating nuclear 
effects. 

We sum the different error estimations linearly, since there might be cor- 
relations which we are not aware of. This assumption seems wrong mainly for 
the radiative corrections, which seem independent of the other contributions. 
Had we took a statistical sum for those, the error would have decreased to 
about 0.7%. However, a conservative estimation is appropriate due to the 
importance of the conclusions. 



6. Discussion 



The conservative error estimation still allows rather interesting conclu- 
sions. First, one notices that the calculated capture rate agrees with the 
experimental measurement T(fi~ + 3 He — > + 3 H)^ = 1496(4) Hz. We 
thus conclude that the hybrid approach used in the current work, which is 



usually named EFT* [22|, accurately predicts the capture rate. 

From this, one practical conclusion can be derived concerning the ability 
of the method of calculation to predict weak reaction rates in astrophysical 
environment, in particular in supernova environment. These reactions are 
usually unreachable experimentally, and their calculation includes momen- 
tum transfer of few tens of MeV. The currents used in these calculations 
are based on extrapolation of f3 decay surveys and theoretical consideration. 
By using the same currents to predict muon capture rates, the extrapola- 
tion becomes interpolation. This conclusion is also deduced by Zinncr et. 



a/.[32j, who used RPA approach to calculate total muon capture rates in 
heavy nuclei. 

However, the most interesting result concerns the weak form factors of the 
nucleon. In order to constrain the induced pseudoscalar and second class form 
factors, we take the following approach. In each case, we set all the other 
form factors to their nominal value, and change this form factor in a way 
which keeps an overlap between the experimental rate and the theoretically 
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allowed rate. The nominal value of the form factor is set to reproduce the 
experimental measurement. 

The resulting constraint on the induced pseudoscalar form factor is: 

gp (q 2 = -0.954m;;) = 8.13 ±0.6, (9) 

in very good agreement with the HB^PT prediction of Eq. ([6]). Together 
with the MuCap results Q, gp{q 2 = -0.88m^) = 7.3 ± 1.2, this is a great 
success to the HB^PT prediction. 

A second conclusion concerns the contribution of second class currents. 
The axial G-parity breaking term was predicted, based on QCD sum-rules to 
be = —0.0152(53) [30] - Using this prediction does not change the result 
of the current calculation significantly (about 0.15%). Our constraint has a 
much larger error bar than this calculation, and agrees with a vanishing form 
factor: 

— = -0.1 ±0.68. (10) 
9 A 

Wilkinson [3l| has collected the experimental data to get \g t \ < 0.3 at 90% 



CL, which is a factor of 2 better than the current limit. One has to still 
consider that the G-parity breaking terms can also excite mesonic currents 
in the nucleus, which were not taken into account in the current discussion. 

To date, the tightest constraints on the vector G-parity breaking term, 
6J}[, related with the Conserved Vector Current (CVC) hypothesis, were 
made using a survey of superallowed + — ^ + /3 decays. This leads to a 
value of g s = 0.01 ± 0.27 [4j. Our calculation puts a much tighter limit, 
which can be considered experimental, on this form-factor: 

9s = -0.005 ±0.040,. (11) 

Thus, the current limit provides the tightest constrain on g s , and still agrees 
with CVC. 

Summarizing, we have calculated the rate of the weak process 3 He(/i~, z/ M ) 3 H. 
The calculation predicts a capture rate of V — 1499 ± 16 Hz, in accord with 
the measured rate T = 1496 ± 4 Hz. The error estimation has two main 
sources, uncertainties in the experimental triton half life, and in the calcu- 
lation of radiative corrections to the process. The nuclear wave functions 
are calculated, using the EIHH method, with the phenomenological nuclear 
forces AV18±UIX. The hadronic currents within the nucleus are derived from 
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HB%PT. Their low energy constants are constrained from low energy pion- 
nucleon scattering, and from the triton half-life. As a result, the calculation 
has no free-parameters. The induced pseudoscalar form factor is constrained 
to ±8%, and agrees with HB^PT prediction [24]. We show that this predic- 
tion is consistent with vanishing second class terms. The CVC hypothesis is 
confirmed to a new limit \g s \ < 0.045. The calculation shows that nuclear ab 
initio calculations of muon capture process can have percentage level accu- 
racy, and thus can be used as a quantitative test for the weak structure of 
the nucleon and other properties of QCD at low energy. 
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